Orbits in a central force field: Bounded orbits 
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The nature of boundedness of orbits of a particle moving in a central force field is investigated. 
General conditions for circular orbits and their stability are discussed. In a bounded central field 
orbit, a particle moves clockwise or anticlockwise, depending on its angular momentum, and at the 
same time oscillates between a minimum and a maximum radial distance, defining an inner and an 
outer annulus. There are generic orbits suggested in popular texts displaying the general features of 
a central orbit. In this work it is demonstrated that some of these orbits, seemingly possible at the 
first glance, are not compatible with a central force field. For power law forces, the general nature 
of boundedness and geometric shape of orbits are investigated. 
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I. INTRODUCTION 

The central force motion is one of the oldest and widely 
studied problems in classical mechanics. Several familiar 
force-laws in nature, e.g., Newton's law of gravitation, 
Coulomb's law, van-der Waals force, Yukawa interaction, 
and Hooke's law are all examples of central forces. The 
central force problem gives an opportunity to test one's 
understanding of the Lagrange's equation, Hamilton's 
equation, Hamilton Jacobi method, and classical pertur- 
bation. It also serves as an introduction to the concept of 
integrals of motion and conservation laws. We need only 
to appeal to the principles of conservation of energy and 
angular momentum to describe the nature and geometry 
of the possible trajectories in central force motion. Most 
books in classical mechanics^S* 3 ^*^, treatise^, and ad- 
vanced texts&iii discuss the central force problem. In this 
article we present some interesting features of bounded 
orbits in a central field. 



A. Kepler's Laws 

One of the most remarkable discoveries in the history of 
physics is that of Keplerian orbits. A tremendous wealth 
of data on planetary positions was collected by Tycho 
Brahe and Johannes Kepler after detailed observation 
spread over several decades. After a thorough analysis 
of this data Johannes Kepler formulated three empirical 
laws that described and correlated the motion of the five 
planets then known: 

1. Each planet moves in an elliptical orbit, with the 
sun at one of its foci. 

2. The radius vector from the sun to each planet 
sweeps out equal areas in equal times. 

3. The square of the periods (T 2 ) of the planets are 
proportional to the cube of the lengths of the cor- 
responding semimajor axes (a 3 ). 



B. Newtonian Synthesis 

Almost 100 years later Newton realized that the plan- 
ets go about in their nearly circular orbits around the sun 
under the influence of the same force that causes an ap- 
ple to fall to the ground, i.e., gravitation. Newton's law 
of gravitation gave a theoretical basis to Kepler's laws. 
Kepler's laws can be derived from Newton's law of grav- 
itation; this is often referred to as the Newtonian syn- 
thesis. Kepler's first and third laws are valid only in the 
specific case of inverse square force. There are, however, 
certain general features which are observed in all cen- 
tral field problems. They include (i) certain conserved 
quantities (energy, and angular momentum), (ii) planer 
nature of orbits, and (iii) constancy of areal velocity (Ke- 
pler's second law). A large class of central forces allows 
circular orbits (stable or unstable), bounded orbits, and 
even closed and periodic orbits. Certain common charac- 
teristics about the generic shapes of bounded orbits can 
also be ascertained. 



II. EQUATIONS OF MOTION AND THEIR 
FIRST INTEGRALS 

A. Central field orbits: confinement in a plane 

The central force motion between two bodies about 
their center of mass can be reduced to an equivalent one 
body problem in terms of their reduced mass m and their 
relative radial distance r. Hence in this reduced system, 
a body having the reduced mass moves about a fixed 
center of force. 

Consider the motion of a body under a central force, 
F = F(r) = f(r)r with the origin as its force center. The 
potential V(r) from which this force is derived is also a 
function of r alone, F = — VV, V = V(r). 

On account of the central nature of the force, the me- 
chanical properties of the body do not vary under ro- 
tation in any manner around the center of force. Let 
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the body be rotated through an infinitesimal angle 56, 
where the magnitude 59 is the angle of rotation while 
the direction is that of the axis of rotation h. The 
change in the radius vector from the origin to the body 
is | 8r \— rsin(9)59, with i5r being perpendicular to r 
and 59. Hence 5r = 56 x r. The change in velocity is 
similarly given by 5v = 56 x v. The Lagrangian of the 
system is a function of r and r. The motion is governed 
by the Lagrange's equation, 

d fdL\ dL 
dt \ dr J dr 







If we now require that the Lagrangian of the system re- 
main invariant under this rotation we obtain 

5L = — ■ 5r + — ■ 5r = . 
or or 

One can define generalized momentum as, 

dL 

p = ^ 

From the Lagrange's equation we get, 

P ~di \jjr~) ~ ~dr~ 

Replacing dL/dr by p and dL/dr by p in the equation 
for 5L we get, 

p-<50xr + p-<50xr = O. 
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Area = 1/2 (r) ( r d G ) 



d0 



FIG. 1: Area swept by radius vector 



C. First integrals and conservation laws 

The canonical momentum corresponding to 9 is called 
the angular momentum (or rather the magnitude of the 
angular momentum that we discussed before), 

Pe = — r = mr 9 = I 

de 

As 9 is a cyclic coordinate, i.e., the Lagrangian is inde- 
pendent of 8, this angular momentum is conserved. This 
can be shown from the Lagrange's equation for 9. 







As 56 is arbitrary we conclude 1 = r x p is a conserved 
quantity. 1 is called the angular momentum of the sys- 
tem. Since 1 is a constant and is perpendicular to r it 
follows that the radius vector of the particle lies in a 
plane perpendicular to 1. This implies that the motion 
of the particle in a central field is confined to a plane. 



d_ fdL\ _ dL _ 
dt Vdl) ~ ~d9 



Pe 



dt 



{mr 2 9) = 0. 



The corresponding integral of motion is, 
mr 2 9 = I 



(4) 



B. Lagrangian and equations of motion 

As the motion in a central force field is confined to a 
plane, it suffices to use plane polar coordinates. One may 
write the Lagrangian of the particle as, 



T 1 -2 

L = -mr 
2 
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V(r) = ^m(f 2 +r 2 9 2 ) -V(r) . (1) 



We assume the center of force to be at the origin. The 
coordinates of the body of mass m undergoing the central 
field motion are given by (r, 9) . 

The Lagrange's equation for the 9 and r coordinates 
are given respectively by, 





fdL\ 


dL 


-1 

dt ' 


KdOj 


~ ~d9 




(dL\ 


dL 


-1 

dt 1 


V dr 


dr 







(2) 
(3) 



and this / can easily be shown to be the magnitude of the 
angular momentum vector 1 = r x p. 

This conservation law is essentially equivalent to Ke- 
pler's 2nd Law: 

The elementary triangular area swept by the radius 
vector in an infinitesimal time interval dt is, 

dA = ^r(rdff) 

hence it follows from Eq. (4) that the rate of areal sweep 
is a constant. 



dA 



1 



dt 2 r 2m 



(5) 



There is another first integral of motion associated 
with the Lagrange's equation for the r coordinate. 

d . . . - 2 dV n 

— (mr) — mr8 + — — = 
dt dr 
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FIG. 2: Constancy of angular momentum in an inverse square 
potential. Total energy is 60% of the minimum of effective 
potential V(r). 



FIG. 3: Constancy of total energy in an inverse square poten- 
tial. Total energy is 60% of the minimum of effective potential 
V(r). 



The force in terms of the potential (conservative force) 
is given by, f(r) = —dV/dr, whence the above equation 
becomes 



mr — mrb 



Using the first integral of motion, one can convert this 
second equation into an equation for r alone. 

d / 1 I 2 
mr = -— V+- j 

dr \ 2 mr z 



Integrating we get 



1 



1 I 2 



-mr 



2 mr' 



+ V(r) = E 



(6) 



where E is a constant of integration, called the energy. 
This is the law of conservation of total mechanical energy 

It is interesting to note that for motion in a general 
force field, I = mr 2 9 remains invariant even though r and 
9 vary with time (see Fig. 2). Similarly E — mr 2 /2 + 
I /(2mr ) + V(r) remains a constant even though f and r 
and hence mr 2 /2 and V(r) + I 2 / (2mr 2 ) each varies with 
time (see Fig. 3). 

Lagrange's equations are two second order Ordinary 
Differential Equations (ODE) in r and 9. However they 
decouple, i.e., each equation is expressible in terms of ei- 
ther r or 9. On integrating each equation once we get 
the first integrals of motion namely, the total mechanical 
energy and the angular momentum. A further integra- 
tion will yield the complete solution to the problem. This 
second integration introduces two more constants of in- 
tegration namely, the initial radial (r ) and angular (9 ) 
positions. 



D. Equation for the orbit 

From the equation giving energy as a first integral of 
motion, we get the expression for radial velocity, 




On integration we get, 



t = 



E-V 



dr 



I 2 



2mr 2 



(7) 



(8) 



, ro y/2/m(E - V - I 2 /(2mr 2 ))' 

This relation can be inverted to give r as a function of t, 
r = r(t). The other first integral gives, 

• I 



mr 



Substituting the expression for r as a function of t, and 
integrating, 



9-9 Q = 



I 



dt 

W 



(9) 



These two expressions for r(t) and 6(t) express the equa- 
tion of the orbit for the particle in a central field in terms 
of time t as a parameter. 

Instead of expressing the orbit parametrically in terms 
of t, one often wants to express the orbit directly as an 
equation connecting r and 9. Such an equation may be 
obtained by eliminating t from the above expressions for 
f and 9. 



d6 
dr 



l/(mr 2 



sj2jm(E - V(r) - l 2 /(2mr 2 )) 
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On integration this yields, 
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l/(mr 2 ) 



y/2 /m(E - V(r) - l 2 /(2mr 2 )) 



dr 



For understanding the qualitative nature of motion in 
a central field one looks at the equivalent one dimensional 
problem. 



2 I 2 
m v w 2mr 2 ' 



-(E-V(r)). 



We call V the effective potential, introduced to make 
the problem similar to that of a particle moving in a 
one dimensional potential field. The effective radial force 
(f(rj) is connected to the effective radial potential (V(r)) 
by the expected relation, 



/» = - 



dV(r) 
dr 



(10) 



At a point where the effective potential V(r) equals the 
energy E, the radial velocity vanishes (f = 0). In one 
dimensional motion this corresponds to a particle com- 
ing momentarily to rest, and having zero kinetic energy. 
However in the case of central field, the motion is not 
really one dimensional, and even for r = 0, the particle 
is not at rest (v = r60), and it has a non-zero kinetic 
energy ((l/2m)r 2 8 2 ). 

For the inverse square force, as in the case of gravita- 
tion, we have f(r) = —k/r 2 and V(r) = —k/r. Effective 
potential (see Fig. 4) is given by, 



V(r) 



2mr 2 



(11) 



The following properties of the effective potential are 
easily noted: 

1. V(r) = l/r 2 {l 2 /2m - kr), as r -> the term 
within bracket is essentially l 2 /{2m) and hence 

lim r _^ V( r ) +°° 

2. linv^oo V(r) = 

3. V(r) = —l/r(k — l 2 /2mr 2 ) and for large values of 
r, l 2 /2mr 2 is negligible compared to k, hence V(r) 
has a negative value. 

4. At r* = l 2 /2mk the function V(r) intersects the r 
axis, i.e., V(r*) = 0. 

5. V(r) reaches a minimum at ro = 2 ■ r* = l 2 /mk. 
With dV/dr\ ro = and d 2 V/d 2 r\ ro > 0. 

For total energy E = V(r ) the particle has zero radial 
velocity at r — r , and no other radial position is phys- 
ically accessible, since r becomes imaginary at r ^ ro. 
This corresponds to circular motion. 



(a) circular orbit 

(b) elliptic orbit 

(c) elliptic orbit of higher eccentricity 
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FIG. 4: Effective potential for an inverse square force 

For > E > V(ro) there is a range of radial positions 
(r ma x > r > r min ) for which E > V(r). The particle can 
move in this range of r with varying r. The radial velocity 



vanishes at the end points r max and 



where the 



energy E equals the effective potential V and r reaches a 
maximum at r — r The points r max and r min are called 
the turning points. The central field particle cannot move 
beyond these points, as the energy becomes less than the 
effective potential, and the expression for radial velocity 
(f) turns imaginary. 

E < V(ro) is a physically impossible situation, since 
no (radial) position is physically allowed for the particle. 

For E > we have an unbounded motion, where the 
particle can fly off to infinity. Thus in the case of inverse 
square force field we can have bounded (E < 0) or un- 
bounded (E > 0) motion depending on the energy of the 
particle. In particular for E — it is parabolic and for 
E > it is hyperbolic. 



III. EXISTENCE AND STABILITY OF 
CIRCULAR ORBITS FOR CENTRAL FORCES 

At all positions other than where the effective potential 
is a minimum or a maximum, we have a net effective force 
f(r) — —dV/dr ^ 0. When the total energy E is not 
equal to the minimum of effective potential V(r), at the 
points of instantaneous zero radial velocity the particle 
is pushed away. If the effective potential has a minimum 
and energy is greater than that minimum then the radial 
distance has a lower and an upper bound and the particle 
moves between these radial limits. On the other hand if 
the effective potential has a maximum, and the energy is 
less than that maximum, the effective force pushes the 
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FIG. 5: Spiralling orbits, f(r) ~ -1/r 



particle away from the positions of zero radial velocity 
in an inward or an outward spiral, ff the total energy 
is equal to the maximum or minimum of the effective 
potential, the system can stay with zero radial velocity, 
and hence move in a circular orbit. 
The condition for circular orbit is, 



dV 
dr 



dV 
dr 



I 2 



whence we get, 



dV 



(12) 



(13) 



The negative sign on the right hand side clearly shows 
that the force must be attractive. The particle moves in 
a circular orbit, since the force of attraction due to the 
central field provides the necessary centripetal force. 

For a given attractive central force ( f(r) and V(r) 
given ) it is possible to have a circular orbit of radius 
ro provided the angular momentum and energy of the 
particle are given by, 



I 2 = 
E = 



-mr%f(r ) 



V(r ) 



2mrQ 



(14) 
(15) 



However when the effective potential has a maximum 
the system is in an unstable circular orbit. A small devi- 
ation from this radial position causes the orbit to be un- 
bounded. The effective force that comes into play makes 
the particle move away from the position of zero radial 
velocity in an inward or outward spiral (Fig. 5). When 



the effective potential has a minimum, the effective forces 
cause the particle to remain in a bound orbit confined 
in an annular space. For small deviations the annular 
radii are nearly the same, and hence the orbit is close to 
a circle. One can understand the stability question by 
studying the forces (restoring or unsettling) that come 
into play when the system is moved infinitesimally from 
the position of circular orbit. 
The circular orbit is stable if, 



d 2 v 

dr 2 



> 0. 



(16) 



That is equivalent to, 

d 2 V ( df 3l 2 



dr 2 
— Ir 

dr r " ~~ mr 4 



< 



dr mr 4 
3l 2 



> 



Using Eq. (13), 



8_l 

dr 1 



< 



3/M 



(17) 



IV. NATURE OF BOUNDED ORBITS 

General bounded motion has both lower and upper 
bounds. It means that the particle cannot approach 
nearer than some minimum or move farther than some 
maximum distance. One has to remember that the angu- 
lar velocity has a constant sign, same as that of the con- 
stant angular momentum, throughout the motion. How- 
ever its magnitude decreases with increase in the radial 
distance (~ 1/r 2 ). Together with the angular motion, the 
radial distance changes from r max to a r m j n , then back 
to r max and so on. From this general nature of motion 
two families of generic orbits are suggested in popular 
textsii (Fig. 6, Fig. 7). However it can be shown that 
the generic types shown in the later figure (Fig. 7) are 
not feasible for any attractive potential. 

We consider below a general bounded orbit confined in 
an annular region. We would like to investigate the na- 
ture of the orbit close to the point where it touches the 
inner or the outer annulus. Let us choose the reference 
line (polar line) of the coordinate system such that the 
orbit touches the annular ring at 9 — 0, or more specifi- 
cally at r = R, 9 = 0. Since at this point P(r =11,9 = 0) 
the orbit is at its closest (or farthest) approach from the 
pole, (dr/dt)p = 0. As r is a function of 9, using Eq. (4) 
we get for motion along the trajectory, 

dr I dr 
dt mr 2 d9 



hence 



= 



(18) 



9=0 




-1.5 




FIG. 6: Generic orbits in central field (allowed) 



FIG. 7: Generic orbits in central field (not allowed) 



The equation of motion for r is given by, 

mr — mrO 2 = f(r) 
From which we get r"{9) along the trajectory, 



I d ( I dr 



rn 



mr 2 dO \ mr 2 d6 



mr- 



TflT 



f(r) 



= —f(r) 
dO 2 I 2 JV) 

mr 4 
= r + —.f(r) 



One can expand r{9) in a Taylor series in 9, and remem- 
ber that r'(0) = 0, 



r(0)traj = r(0)+r'(0)6 + r"(0)-+h.o. 

= r(0) + ^/>(0))£ + A. O . 

= r(0) + (r(0) + l^/(r(0))) ^ + h.o. 



Consider a tangent to the annulus at the point P, Fig. 
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tangent 



trajectory 
(allowec 




tangent 



trajectory 
(not allowed) 



FIG. 8: Curvature of the trajectory in comparison to the 
tangent to the inner annulus 



FIG. 9: Curvature of the trajectory in comparison to the outer 
annulus 



r(0)tan=r(0)+r(0)- + h.o. 

Notice that the constant or the 9 independent terms in 
r{0)traj and r(9) tan are equal, and the first leading power 
of 9 is 9 2 in both the cases. The coefficient of 9 2 (cur- 
vature) will determine the nature of the trajectory in 
reference to the tangent. 

We first study the curvature of the orbit in reference 
to the tangent. The force being attractive /(r(0)) < 0, 
for small angular distance (A(0)) away from the point P, 



we study the orbit near the point P where it touches the 
outer annulus. 



which can be expanded to give r{9) tan along the tangent, hence 



r(A6) traj < r(0) = r„ 



(20) 



This confirms that the outer annulus is indeed the outer 
bound of the trajectory. The nature of the orbit near the 
point of contact P at the outer annulus is shown in figure 
9. 



V. BOUNDED ORBITS FOR THE POWER LAW 
CENTRAL FORCE 

A. Existence of stable circular orbit 



I (I 



(19) 



which means that the orbit bends more sharply than the 
tangent and stays closer to the inner annulus for nearby 
points (Fig. 8). Hence the second type of orbits shown 
in some texts (Fig. 7) are not possible for central forces. 
For the outer annulus the analysis with respect to tan- 
gent is not very meaningful as one can set up an even 
stronger bound for its orbit, directly considering the po- 
tential V(r). In this case the effective potential has a 
positive slope with respect to the radius, and hence a 
negative effective force (/(r) < 0). The orbit should not 
only remain nearer than the tangent, but even nearer 
than the outer annulus. This condition is confirmed if 



For the case of a power law central potential 



Sir) = 



nk 

r n+l 



From the stability condition Eq. (17), 



we find, 



Or 1 



(n + l)nk 



< - 



3/(r ) 



„n+2 



< 



-3 — 



.71+ 1 



n < 2 



1 

ro 



(21) 



(22) 
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V(r) 
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FIG. 10: Effective potential for f(r) ~ -1/r 4 force 

Hence the circular orbit is stable for an attractive power 
law potential that varies slower than inverse square (or a 
force that varies slower than inverse cube). 



For the point of extremum 

„ 2 amn 

Finding the second derivative of V with respect to 
r at the point r , 

d 2 V anjn+V) 3l 2 

dr 2 r n+2 mr 4 

1 / an(n+l) 3l 2 \ 

T 4 \ r "-2 m J 

If, A 2 3Z 2 \ 

= I -(« + !)- + — 
Tq \ mm/ 

1 I 2 

= --(2-n)<0 

Hence it is a point of maximum. For E > V(r^) we 
always have an unbounded orbit. For E < V(tq) the 
orbit is semiboundcd, bounded above or bounded below, 
according to its initial state. The particle either spirals in 
or spirals out. For E = V(ro) we get an unstable circular 
orbit. 



2. V(r) = -a/r n where n = 2 



B. Study of boundedness of orbits 

1. V(r) = -a/r n where n > 2 

Consider the case when n > 2. We have the effective 
potential, 



V(r) = - 



(23) 



1. V(r) = -l/r n (a - I 2 ■ r"- 2 /(2m)), and as r -» 
we may neglect I 2 ■ r n ~ 2 /(2m) in comparison to a, 
thus lim r ^o V( r ) = —oo- 

2. linv^oo V(r) = 0. 

3. For large but finite r, V(r) = -l/r 2 (a/r n ~ 2 - 
I 2 /(2m)), and a/r n ~ 2 is negligible compared to 
I 2 /(2m), and hence V(r) is positive. 

4. V(r) intersects the r axis at r* where ( r *)"~ 2 = 
2am/; 2 . 

5. V(r) has a maximum at r where r^ 2 = anm/l 2 . 

For a maximum point of V we should have 
dV/dr \ ro = and d 2 V/dr 2 \ ra < 0. 



dV an I 2 



The effective radial potential is V(r) = —a/r 2 + 
l 2 /(2mr 2 ). This is essentially an attractive or repulsive 
inverse square term. 

a I 2 
V(r) = - — , a > if a > — 



a < if a < 



2m 
2m 



dr r n+l mr 3 



= 



This potential cannot give circular orbit ever. If the ef- 
fective potential is attractive and E < it has an upper 
bound of radial distance. A typical orbit would therefore 
be an inward spiral. For a repulsive effective potential 
we will have outward spiral moving to infinite radial dis- 
tance. 



3. V(r) = -a/r n where 2 > n > 

In this case the effective potential is, 

a I 2 

V(r) + ^ (24) 

v > r n 2mr 2 v ; 

We have the following properties of V(r) 

1. V(r) = l/r 2 (l 2 /2m-a/r n ~ 2 ) and as r — > we may 
neglect a/r" -2 in comparison to I 2 /2m, and thus 
lim r _^ V( r ) +oo. 

2. lim^oo V(r) = 0. 
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FIG. 11: Effective potential for f(r) ~ -1/r 2 force 
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FIG. 12: Effective potential for f(r) ~ — 1/r force 



3. V(r) = -l/r n (a - Z 2 /2mr 2 -") and for large but 
finite values of r, l 2 /2mr 2 ~ n is negligible compared 
to a, hence V'(r) is negative. 

4. V'(r) intersects the r axis at r* where (r*) 2 ~™ = 
l 2 /(2am). 

5. V'(r) reaches a minimum at r where r 2 ~™ = 
Z 2 /amn. 

For a minimum point of V we should have 
dV/dr \ ro = and d 2 V/dr 2 \ ra > 0. 



<9F 
dr 



an 



I 2 



= 



For the point of extremum 

I 2 



ro 



2-n 



amn 



Finding the second derivative of V with respect to 
r at the point r , 



d 2 V 
dr 2 



an(n+l) 3l 2 



„n+2 



an(n+l) 3/ : 



, / ^ l Zl 2 



m 



rn 



1 I 2 

= --(2-n) >0 



to 



Hence it is a point of minimum. For this case the orbit 
can be bounded or unbounded depending on the total 
energy E of the particle. 



4- V(r) — a\nr 

This corresponds to one over r force f(r) 
effective potential is given by, 



V(r) = ah\r + 



l 2 



-a/r. The 
(25) 



2mr 2 

We have the following properties of V(r) 

1. lim r ^o V(r) — > +oo. 

2. lim r ^oo V(r) — > +oo. 

3. There is no point of intersection with the r axis, 
and V(r) is always positive. 

4. V(r) reaches a minimum at ro = I / \/am. 

For a minimum point of V we should have 
dV/dr \ ro = and d 2 V/dr 2 | ro > 0. 



dV 
dr 



a 
r 



I 2 



m,r" 



for the point of extremum 
I 2 

2 1 

r o = 

am 

The second derivative of V" at r 
d 2 y a 3/ 2 



dr 2 



4 I 



3Z 2 



= ^ • 2a > 
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V(r) ~ r"2 



V(r) 



l A 2/(2m r A 2) 
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FIG. 13: Effective potential for f(r) ~ — r force 



The second derivative of V with respect to r, 

„n-2 

1 4 

mr 

3F 



d2 V , n -2 3/ 2 

— - an(n-l)r n 2 + 

ar z mr* 



1 



an(n - l)r"+ 2 + 



rn 



i ; 2 

(2 + „) >0 

r 4 m 



Hence it is a point of minimum. The orbit is always 
bounded. 

The findings for the general power law potential can 
be summarized as follows, 



and 



V(r) — sign(n)ar n 
f(r) = -absi^ar 71 ' 1 



V(r) = folnr 

m = 4 



(27) 
(28) 

(29) 
(30) 



Hence it is a point of minimum. The orbit is always 
bounded. 



5. V(r) — ar n where n > 

For this potential the effective potential is given by. 

I 2 



TABLE I: Dependence of boundedness on power law. 



V(r) = ar n + 



2mr 2 



n>0 



(26) 



We have the following properties of V(r) 

1. lim r ^o V( r ) — * +°°- 

2. limr^oo V{r) — > +oo. 

3. There is no point of intersection with the r axis, 
and V(r) is always positive. 

4. V'(r) has a minimum at r , Tq +2 = l 2 /amn. 

For a minimum point of V we should have 
dV/dr | ro = and <9 2 V7cV 2 | ro > 0. 



= anr 



n-l 



dr mr 
for the point of extremum 

I 2 

r n+2 - — 
amn 



n of V(r) ~ ar n 


nature of boundedness 


)-oo,-2( 


always unbounded 


-2 


spiralling orbit 


)-2,0( 


bounded or unbounded depending on E 


(lnr) 


always bounded 


)0,+oo( 


always bounded 



C. Stable bounded orbit, geometric shape 

The stable bounded orbits in a power law central field 
can have the following forms 

1. V(r) = -a/r n , < n < 2 

2. V(r) = blogr 

3. V(r) =a-r n , n>0 

In all the above cases the derivative of V(r) with respect 
to r is always positive, and hence the force is necessarily 
attractive (/(r) = -dV(r)/dr < 0). 

In the first case the effective one dimensional potential 
V(r) goes to infinity as r 0. V(r) has a minimum at 
some r = ro, and it has a negative slope for all r < tq, 
and positive slope for all r > r . For E = V(r ) we get 
stable circular orbit, and for V(r ) < E < 0. the orbits 
are still bounded. For cases (2) and (3) the circular orbit 
is stable, and orbits are always bounded for any energy. 
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TABLE II: Dependence of stability of circular orbits on power 
law. 



n of V(r) ~ ar n 


stability of circular orbits 


)-oo,-2) 
)-2,oo( 


unstable 
stable 



VI. CONCLUSION 



orbits are sometime mistakenly identified with diagrams 
of the form shown in Fig. 7—. The above expression in 
fact corresponds to a class of orbits that look generically 
like those shown in Fig. 6. 

The generic features of central force orbits discussed 
here have been verified by computer simulation for a large 
class of central force fields. The figures shown here (Fig. 
2, 3, 5 and 6) were generated by these simulations. Stu- 
dents interested in studying and generating such orbits 
will find Ref.12 helpful. 



Existence of bounded orbit for a large class of attrac- 
tive central field has been discussed. The generic nature 
of central field bounded orbits is analytically derived. 
Certain class of these orbits (Fig. 7) presented in popular 
texts-, are shown to be non-feasible. 

Small deviation from circularity in the case of central 
field is often expressed in terms of inverse of radial dis- 
tance (u — 1/r). 

u = uq + a ■ cos(/30) 

r = 

1 + a ■ ro ■ cos([36) 

where ro = 1/uq. It is interesting to note that these 



Acknowledgement 

The authors wish to express their indebtedness to the 
well known texts by Goldstein^, Landau^ and Arnold 9 . 
They also acknowledge their teachers in related graduate 
courses at Stony Brook, Prof. Max Dresden, Prof. A. S. 
Goldhaber and Prof. Leon A. Takhtajan. 

Authors gratefully acknowledge the encouragement re- 
ceived from Prof. Shyamal SenGupta of Presidency Col- 
lege, Calcutta. The material presented here was used in 
a graduate level classical mechanics course at Jadavpur 
University during 1998-2001. SR wishes to thank his stu- 
dents for stimulating discussions. 



* Electronic address: subho@juphys.ern et.in| 

1 H. Goldstein, Classical Mechanics, Addison- Wesley, 1980. 

2 L.D. Landau, E.M. Lifshitz, Mechanics, Pergamon Press, 
1976. 

3 K.R. Symon, Mechanics, Addison- Wesley, 1971. 

4 D.T. Greenwood, Classical Dynamics, Prentice Hall, 1977. 

5 J.L. Synge, and B.A. Griffith, Principles of Mechanics, Mc- 
GrawHill, 1970. 

6 A. Sommerfeld, Mechanics, Academic Press, 1952. 

7 E.T. Whittaker, A Treatise on the Analytical Dynamics of 
Particles and Rigid Bodies, Dover, 1944. 



L.A. Pars, Introduction to Dynamics, Cambridge, 1953. 
9 V.I. Arnold, Mathematical Methods of Classical Mechanics, 
Springer Verlag, 1989. 

10 R. Abraham, and J.E. Marsden, Foundations of Mechan- 
ics, Benjamin-Cummings Publ.Co.Inc, 1978. 

11 H. Goldstein, Fig. 3-7, Fig. 3-13 and Eq. 3-45 in Classical 
Mechanics, Addison- Wesley, 1980. 

12 H. Gould, J. Tobochnik, An Introduction to Computer- 
Simulation Methods, Addison- Wesley, 1988. 



